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Abstract 

Using a liquid-state approach based on Ornstein-Zernike equations, we study 
the behavior of a fluid inside a porous disordered matrix near the liquid-gas 
critical point. The results obtained within various standard approximation 
schemes such as lowest-order 7-ordering and the mean-spherical approxima- 
tion suggest that the critical behavior is closely related to that of the random- 
field Ising model (RFIM). 
Pacs numbers: 05.70.Fh, 61.20.Gy, 64.60.Fr 



1 



Typeset using REVTgX 



The interpretation of experimental studies of the phase behavior of fluids (binary mix- 
tures [0] or ^He 0) adsorbed in disordered porous materials such as silica gels or Vycor 
remains highly controversial. Some results lend support to the view that preferential ad- 
sorption of one component and randomness of the pore network give rise to a random-field 
Ising-like behavior near the critical point, as suggested by Brochard and de Gennes On 
the other hand, it has been argued that the experiments in Vycor, which is a glass with 
a fairly low porosity, can be interpreted in terms of wetting in a confined geometry, with 
no randomness 0. Clear-cut conclusions are somewhat obscured by the metastable and 
non-equilibrium effects which are often present in the experiments. However, there is a need 
for more realistic continuum descriptions which encompass both randomness, confinement, 
and connectivity between pores. The goal of this letter is to present a liquid-state approach 
describing a single component fluid in a rigid disordered matrix and, starting from this mi- 
croscopic theory, to give some insight on the critical properties of the fluid and their relation 
to the random-field Ising model (RFIM) |p. 

Following the formalism recently proposed by Madden and Glandt [^, we consider the 
matrix-fluid system as a special binary mixture in which the matrix is treated as a rigid 
set of obstacles, obtained by quenching an equilibrium configuration of particles (species 0), 
with the fluid particles (species 1) at equilibrium in the presence of (and in interaction with) 
the quenched particles. The structure of this quenched- annealed (Q-A) mixture is described 
by the total pair correlation functions hij{r) which satisfy a set of coupled integral equations, 
the so-called replica Ornstein-Zernike (ROZ) equations [0. As in the usual OZ formalism [§, 
these equations relate the functions hij{r) to the corresponding direct correlation functions 
Cij{r) which are expected to be short-ranged (i.e, J d'^v Cij{r) is finite), even at the critical 
point. The exact form of the ROZ equations can be most easily derived by using the 
replica trick which replaces the study of the original system by that of an equilibrium 
mixture composed of s copies of the fluid interacting with the matrix. It should be stressed 
that, in contrast to the standard replica treatment of quenched disorder, there is no direct 
interaction between the fluid replicas because the average over disorder (i.e., over all possible 
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configurations of tlie matrix particles) is not performed explicitly. These ROZ equations 
represent the starting point of our treatment and their formal solution in Fourier space is 
given by 

Uk) = r^^^ (la) 
= n (1^) 

[1 - pocoo{k)\[l - piCc{k)\ 

where /(k) denotes the Fourier transform of a function /(r) and pi is the average density of 
species i. In deriving these equations, it has been assumed that there is no replica symmetry 
breaking. As in the case of disordered magnetic systems, one has introduced the connected 
part, hc{r), and the blocked or disconnected part, hb{r) = /iii(r) — hc{r), of the fiuid-fiuid 
correlation function they are equal, respectively, to the s = limits of haa{r) — hajsir) 
and hapir), where a and P ^ a denote fluid replicas. Similarly, Cc(r) and Cfe(r) are defined 
as the connected and blocked parts of cii(r). By symmetry, one has /iqi = /^lo and cqi = cio- 
The above ROZ equations differ from those recently derived for a fluid in presence of 
a quenched, random, external potential obeying a Gaussian distribution [|I^: in Eq. (Id), 



the quenched disorder enters explicitly through a term involving the matrix-matrix and 
matrix-fluid direct correlation functions which also appear in Eqs. (la) and (lb), whereas 



in Ref. |10| it is only contained in the inter-replica pair potential. As a result, the present 
formulation is also able to take into account the geometrical constraint on the fluid due to 
the presence of the matrix (excluded volume and/or confinement effects). 

It has been shown @,|ll|] that Cc(0) is related to Xi? the compressibility of the fluid 
inside the matrix (i.e., the derivative of the fluid pressure Pi with respect to the density or, 
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through the Gibbs-Duhem equation 0, the derivative of the fluid chemical potential), via 
the equation 

A;sTpiXi = [l-pic,(0)]-i (2) 

Therefore, the gas-liquid critical point of the fluid inside the matrix, if it exists, is charac- 
terized as usual by both a diverging compressibility and long-range correlations. However, 
a special feature of Q-A systems is that the disconnected part of the fluid-fuid correlation 
function, and as a consequence /iii(r) itself, may be much longer-ranged that the connected 
part. Indeed, if one defines the exponents r] and fj by hc{r) ~ f-d.+2-ri hb{r) ~ ^^-d+A-rj 
when r — i> oo at the critical point, one finds immediately from Eqs. (Ic) and (Id) that 
fj = 21]. This result rests only on the assumption that the matrix is non-critical, i.e., 
1 ~ PoCoo(O) 7^ 0, and that 4(0) and Coi(O) are finite, a condition generally expected for di- 
rect correlation functions. This relation between r/ and t] is for instance satisfied within the 
standard Ornstein-Zernike approximation which assumes that the direct correlation func- 
tions have the same range as the associated pair interactions. For short-ranged interactions, 
this gives r] = rj = 0, but for long-ranged fluid-fluid interactions decaying as r~^'^'^'^^ for 
r — > oo with a < 2, one obtains rJ = 4 — 2cr = 2t]. Note that the same relation f] = 2t] has 
been proposed for the RFIM |jl2| and seems to be confirmed by recent numerical estimates 



IK. 

To be more specific, we now suppose that the fluid-fluid (resp. matrix-fluid) pair potential 
can be divided into a reference part which includes the short-range repulsive part of the 
interaction and a more smoothly varying long-range attractive part Wii(r) (resp. Woi{r)). 
As a first step, we derive the mean-field (MP) results. As in the bulk fluid case, this can be 
done by considering the lowest 7- ordering approximation which is derived by introducing the 
inverse range parameter 7 in the attractive fluid-fluid potential, Wii(r) = 7'^0(7r), and then 
taking the limit 7^0 (see for instance Ref. 0) . From the graphical expansion of the three 
direct correlation functions Coi(r), Cfe(r) and Cc(r) [Q, it is easy to show that Coi{k) = CQ^{k), 
Chik) = c^{k) and cdk) = c^{k) — {kBT)''^(f){0)6kfi, where 6k,o is the Kronecker symbol and 



the superscript R indicates a quantity calculated in a reference system for which wu{r) = 0. 
The fact that Cf, reduces to its reference part illustrates the fact that the fluid replicas are 
not directly coupled (see above). Using Eq. (2), we immediately obtain a van der Waals 
equation of state, Pi = Pi^^ + |pf0(O). From the thermodynamic relations satisfied by the 
Q-A mixture 0, we can conclude that besides 1] = r] = 0, the critical behavior is described 
by the usual MF exponents, u = 1/2, a = 0, P = 1/2, 7 = 1, 5 = 3. As illustrated 



numerically in Ref. |I^ for a specific system, the critical point is displaced from the MF 
critical point of the bulk fluid because the reference system includes the influence of the 



host matrix (note that in Ref. |jT4| the matrix-fluid interaction is also treated in the MF 
approximation) . 

To go further, some more quantitative statement on the correlations over small distances 
must be made. This can be implemented, within the OZ approximation mentioned above, 
by considering approximate closure relations to the ROZ equations like the mean-spherical 
approximation (MSA) |B|. For bulk fluids, the MSA is known to yield the same critical 



exponents as the spherical spin model We now show that the same is true for a Q-A 
system. To avoid unnecessary technical difficulties, we consider the case of a d-dimensional 
lattice gas. We expect that the behavior in the critical region will be identical to that of 
a real continuum fluid. Lattice cells (or sites) can be occupied either by matrix or fluid 
"particles" and the pair potentials contain a " hard-core" part which excludes multiple 
occupancy of a cell (or site). It must be noted that even for the lattice fluid there is no hole- 
particle symmetry because of the interaction with the matrix. Therefore, there is no exact 
correspondence between this microscopic description of a fluid inside a matrix and the usual 



RFIM (see also Ref. [|T6| for a related discussion). More precisely, the standard equivalence 
between the lattice-gas and the spin 1/2 model shows that we are actually dealing with a 
site diluted Ising model where the spins are coupled to correlated external random fields 
which depend on the site disorder variables Application of the MSA closure in the 

matrix-replicas mixture (i.e. ksTca/sir) = — ^^/^(r) for r 7^ ) gives, in the limit s — > 0, 
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coo(r) = - ,^ s , coi(r) = Coi^r.o - /5w^oi(r) 
Cf,(r) = 0, Cc(r) = c°5r,o - /^w^iilr), (3) 

where Cqj^ and c° are functions of temperature and fluid density, which, owing to the structure 
of the ROZ equations, are completely determined by imposing the "core" conditions h^xiv = 
0) = hii{r = 0) = —1. Note that he and hb do not satisfy themselves the core requirement. 
On the other hand, the prescription for cqo is the direct consequence of the core condition for 
the matrix particles, /ioo(r = 0) = — 1 (for simplicity, we assume that there is no attractive 
interaction between matrix particles). 

For clarity, we first consider nearest-neighbor (n.n.) interactions in a hypercubic lattice. 
Generalizing the analysis developped by Stell [|l5l and Theumann |18] for the bulk lattice-gas, 
we write the inverse Fourier transform of x(k) = 1 -l- pihc(k) as 

'"'""'^(■•) - Gir, K^) ^ r ^ (4) 



ksT ' ' ' (27r)'^ ^2 + 2 ^^t,(l - cos k,] 

where k = (fci, k^) and k = [—2d + ksTil — pic[!)/(— piWn)]^/^ is an inverse correlation 
length which vanishes at the critical point. By replacing in the ROZ equations, we find 



(1-po) 

and 



/loi(r) = fG{r, K ) + woiSrfl (5) 



-iPp^Mir) = < + 2^oi/G'(r, - fG'{r, ^) (6) 
Po(l -Po) 

where G"(r, k^) = 9G(r, and / = /^bTcq^ — woiin"^ + 2d). Obviously, hc{r), hb{r) 

and hoi{r) have the same correlation length and rj = r] = 0. If < 4, no finite critical 
temperature exists since this would lead to the unacceptable result that hb{r) and thus /iii(r) 
do not decay to zero when r oo at the critical point. Using the core conditions together 
with Eqs. (4-6) yields an implicit equation relating k to the thermodynamic variables. In 
absence of hole-particle symmetry, the critical point is found by requiring that {dPi/dpi)c = 
{d'^Pi/dpDc = on the critical isotherm, which from Eqs. (2) and (4) is equivalent to 
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Kc = (9k^/9pi)c = 0. For our system, this yields pic = (1 — po)(l + '2yqpo)/2{l + qpo) with 

q = — 1 — G"(0, 0)/G^{0, 0) > and y = Wqi/wu, whereas the shift of the critical temperature 

due to the presence of the matrix is given by Tc/T^^^'' = (1 — Po)[l + %(! — 2/)Q'Po]/[1 + QPo]] 

this ratio being always smaller than 1, one has Tc < T^"^^^ . The behavior of the correlation 

length in the critical regime can now be obtained by expanding the implicit equation for k 

in terms of AT = T — Tc, Api = pi — pic, and n. Using the asymptotic expansion of the 

Green function G{r = 0, k^) [|1^, we finally obtain for > 4 

G(0,0) ,^ , 1 + gpo ,^ .2 u r^ G' {0,0) , 

ksAT + (Api) = ksTc k 

Wu 1 - Po Wu 

+ JTT—^i^ - ^yfi^'^' + Bn"-' + CK'''HnK + ...] (7) 
4(l + gpo)2' 

where A, B are nonzero constants and C = when and only when d is odd. From this we can 
read off the exponents of the correlation length along the critical isochore (z/) and the critical 
isotherm (e, in the notation of Ref. ||15[), and then, with the help of the compressibility 
equation (2) and the Gibbs-Duhem relation 0, the values of 7 and 6. For d > 6, one 
recovers the MF exponents whereas for 4 < d < 6 one has u = l/{d — A), e = 2/{d — A), 
7 = 2/{d — A) and 6 = d/{d — 4). On the other hand, when y = 1/2, one finds the MF 
exponents as soon as c? > 4. In this latter case, the system is easily shown to be isomorphic 
to the site diluted spin model with no random fields. 

Using the expression of the configurational internal energy, Ui = l/2pit/)ii(0)+poti)io(0) + 
1/2 ^j._^q[Pi 1^11 (r) /ill (r) + 2po?i'io(r)/'-io(r)], and deriving with respect to temperature, one 
finds that the critical exponent a of the fluid specific heat Ci^ is always zero whereas the 
exponent characterizing the singular part of Ciy, if present, is always negative and given 
by = ((i — 6)/((i — 4) for 4 < < 6 and = (6 — d)/2 for d > 6. Finally, the same 



arguments as for the bulk lattice-gas ||T5[ lead to (3 = 1/2 for the coexistence curve exponent. 
These calculations can be easily generalized to non n.n. interactions and to cases where 
woi{r) and Wii(r) do not have the same range. The same exponents are recovered except 
that when tfii(r) dies off like r^^'^^'^^ with cr < 2 one finds u = l/(d — 2a), 'j = a/{d — 2a), 
6 = d/{d — 2a) and as = {d — 3a) /{d — 2a) for 2a < d < 3cr, whereas z/ = 1/(T, 7 = 1,5 = 3 



and tts = (3o" — d)/a for d > Scr. 

Except for the special case which reduces to the site diluted spin model without random 



field, all these exponents coincide with those of the spherical model in a random field []T9 
in particular, we find that the disorder induces a dimensional shift d ^ d — 2 for short-range 
interactions and d ^ d — a ioi long-range interactions. These results are valid even in the 
absence of attractive matrix- fluid interaction = in the case of n.n interactions). The 
fact that the Q-A system and the random-field spin model have the same behavior in the 
critical region does not garantee, however, that their thermodynamics coincide, and a more 
complete study of the MSA phase diagram will be presented elsewhere [|17 . 



To summarize, we have related the critical behavior of a fluid in a porous matrix to 
that of the RFIM by using common approximations of liquid-state theory [|20|. The present 
treatment based on the ROZ equations takes into account both disorder and excluded volume 
effects. The fact that the connection with the RFIM holds also in the absence of matrix-fluid 
attractive interactions suggests that preferential adsorption may not be a crucial ingredient 
for observing a random-field like behavior. Of course, the identity of the critical exponents 
within various approximation schemes is not a proof that the two systems belong to the 
same universality class |^. To go further, one must implement the techniques of the 
renormalization group, which can be considered in connection with liquid-state theory p2 



One must also investigate a possible replica symmetry breaking (RSB) mechanism which 
seems to occur in the RFIM (however, this is not expected at the MF or spherical- 



model levels For that purpose, it is necessary to use a non-linear closure to the ROZ 



equations, such as the hypervertex approximations discussed in Ref. [^. If RSB indeed 
occurs in Q-A systems, the assumption that the direct correlation functions remain short- 
ranged, even in the critical region, may lead to non trivial predictions. 
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